Quantum Computing with Quantum Dots on Quantum Linear Supports 



O 

o 

>>■ 



K. R. Brown, D. A. Lidar,* and K. B. Whaley 
Department of Chemistry, University of California, Berkeley 

Motivated by the recently demonstrated ability to attach quantum dots to polymers at well- 
defined locations, we propose a condensed phase analog of the ion trap quantum computer: a 
scheme for quantum computation using chemically assembled semiconductor nanocrystals attached 
to a linear support. The linear support is either a molecular string (e.g., DNA) or a nanoscale rod. 
The phonon modes of the linear support are used as a quantum information bus between the dots. 
Our scheme offers greater flexibility in optimizing material parameters than the ion trap method, 
but has additional complications. We discuss the relevant physical parameters, provide a detailed 
feasibility study, and suggest materials for which quantum computation may be possible with this 
approach. We find that Si is a potentially promising quantum dot material, already allowing a 5-10 
qubits quantum computer to operate with an error threshold of 10~^. 



I. INTRODUCTION 



The tremendous excitement following the discovery of fast quantum algorithms has led to a proliferation of 
fSJ , quantum computer proposals, some of which have already been realized in a rudimentary fashion. A representative 
CN ■ list includes nuclear spins in liquids and solids Q , trapped ions , atoms in microwave cavities Q , atoms 

in optical lattices jl^], atoms in a photonic band gap material |l2 Jl5] , quantum dots |16-24|, donor atoms in silicon 
p5| , p6t and silicon-germanium arrays Josephson junctions |28p32|, electrons floating on helium | |3^ , electrons 
[ transported in quantum wires p^ , ^ , quantum optics | ]36| , ^ , quantum Hall systems ||3^, and anyons~ p9| , ^ . For 
. critical reviews of some of these proposals see p]-^. To date, no single system has emerged as a clear leading 
' candidate. Each proposal has its relative merits and flaws with respect to the goal of finding a system which is both 
: scalable and fault tolerant Q], and is at the same time technically feasible. In this paper we examine the possibility 
of making a solid state analog of a scheme originally proposed for the gas phase, namely trapped ions. One purpose 
of making such a study is to undertake a critical assessment of both the benefits and the disadvantages which arise 
on translation of an architecture designed for atomic states coupled by phonons, to the corresponding architecture for 
condensed phase qubits. Our proposal uses quantum dots (semiconductor nanocrystals) and quantum linear supports 
(polymers or nanorods) in an ultracold environment. It relies on recent advances in the ability to chemically attach 
^ . nanocrystals to polymers in precisely defined locations. Quantum dots are coupled through quantized vibrations of 
d ' the linear support that are induced by off-resonance laser pulses and information is stored in exciton states of the dots. 
^ Internal operations on exciton states are accomplished using Raman transitions. We provide here a detailed analysis 
that allows evaluation of the merits and demerits of a condensed phase rather than gas phase implementation. 

Semiconductor nanostructures are known as "quantum dots" (QDs) when their size is of the order of or less than the 
bulk-exciton Bohr-radius. In such "zero-dimensional" QDs the electron-hole pairs are confined in all three dimensions 
and the translational symmetry that holds for bulk semiconductors is totally lost. As a result of this quantum 
' confinement the energy-level continuum of the bulk material changes into a discrete level structure. This structure is 
very sensitively dependent on the QD radius and shape, crystal symmetry, relative dielectric constant (compared to 
the surrounding medium), surface effects, and defects. This sensitivity can be used to create and control a wide range 
of optical effects |Q. In general, the term "quantum dot" is used to refer to both "O-dimensional" semiconductor 
structures embedded within or grown on a larger lattice, i.e., lattice bound, and to individual, chemically assembled 
semiconductor nanocrystals |4£]. QDs can be created in a larger crystal structure by confining a two-dimensional 
electron gas with electrodes |47], or by making interface fluctuations in quantum well s [^Sf . A number of promising 
proposals for quantum computation have been made using the lattice-bound dots [p^-^2[. We consider here instead 
the chemically assembled semiconductor nanocrystals. In the remainder of this paper the term QD will therefore be 
implicitly understood to refer specifically to chemically assembled nanocrystals. 

A large amount of theoretical and experimental information about nanocrystal QDs exists. Nanocrystals have 
been studied for their photoluminescence properties, linear absorption properties, and non-linear spectroscopy, using 
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a variety of models and techniques [|49|-p5[. For reviews see, e.g., [p6|-p8[ . These studies clarified the roles of size- 
dependence, lattice structure, surface effects and environment on the exciton spectrum. However, little attention has 
been paid so far to the possibility of using nanocrystal QDs for quantum computing. One reason may be the difficulty 
of coupling nanocrystals. Direct interactions between separate dots are small and difficult to engineer, so that the 
route to scalability is not obvious. In the only other study to date that proposed to use nanocrystals for quantum 
computing, Brun and Wang considered a model of nanocrystals attached to a high-Q microsphere and showed that the 
interaction between QDs can be achieved by using whispering gallery modes of the microsphere to entangle individual 
qubits |2^. One problem with realization of this model is that only a few QDs can be placed on each microsphere. 
Therefore, scalability would depend on the ability to connect the microspheres by optical wires. 

An exciting route to bypass the coupling problem for quantum dots is suggested by the recently demonstrated 
ability to attach QDs to polymers by chemical methods at well-defined locations We show below that at 

sufficiently low temperatures, the QD-polymer system has quantized vibrational modes that can be used to couple 
electronic excitations in quantum dots in a controlled and coherent manner. This "quantum information bus" concept 
derives from the ion trap implementation of quantum computation proposed by Cirac and ZoUer [Q. Ion trap 
schemes take advantage of addressable multilevel ions that are trapped in harmonic wells. The ions are then coupled 
through interaction with their collective vibrational modes This scheme can be extended to any system of 

multilevel quantum objects bound by coupled quantum harmonic oscillators. We apply this approach here to a series 
of nanocrystal QDs attached to a linear support. The excitonic states of the QD act as carriers of quantum information 
which are coupled to the vibrational states of the linear support. A linear support is a one-dimensional material {e.g., 
a stretched polymer or a clamped nanoscale rod) that is connected at each end to a wall. The support is contained 
in either a vacuum or a non-interacting condensed phase matrix such as liquid helium. 

The main advantage of using quantum dots rather than ions is the ability to control the optical properties of 
quantum dots by varying the size, shape and composition of the dot. On the other hand, a disadvantage is that the 
analysis for quantum dots is complicated by the fact that they are complex composite objects and are not naturally 
"clean". For example, defects and surface effects can influence the electronic properties |Q. Our model presupposes 
that nanocrystals which are sufficiently "clean" will ultimately be available, so this puts some severe demands on the 
experimentalist . 

Section || gives an outline of the proposal, describing the basic physics and the formal similarities with the ion trap 



scheme. Section HI describes the physics of the qubits, namely the electronic states of quantum dots, and the quantum 



linear support which provides the information bus between qubits. A summary of the necessary requirements of the 



qubit states is given here. In Section IV we then show how one-qubit and two-qubit operations can be performed in 
this system of coupled quantum dots. Section ^ discusses the feasibility of undertaking quantum logic, with a detailed 
analysis of the constraints imposed by decoherence and physical parameters. Quantitative estimates are made for 
several specific candidate systems in Section followed by conclusions and discussion in Section VIl . 



II. THEORETICAL OVERVIEW 

We outline here the basic elements of the quantum dot-quantum linear support scheme for quantum computation. 
The proposed system consists of semiconductor nanocrystal QDs attached at spacings of several tens of nanometers 
to a quantum linear support (a string or rod). Each QD supports one qubit through a certain choice of excitonic 
states. Single qubit operations are executed by optical transitions between these states. QDs are coupled by the 
linear support in analogy to the ion trap scheme 0|. Thus, one uses detuned laser pulses to excite a phonon of the 
quantum linear support, which can then be used to cause conditional interactions between different dots. The system 
is depicted schematic ally in Fig ure p] . The distance between the quantum dots is assumed to be large relative to their 

|III C| . 



size (see also Section [II B and III Q.) In the presence of external driving fields, the full Hamiltonian can be written 



as the sum of three contributions 

H = Hq + He + Hi , 

where Hq is the free Hamiltonian, He is the coupling Hamiltonian, and Hj is the Hamiltonian describing the inter- 
action between the system and the applied laser fields. 



^ In the original Cirac-ZoUer proposal [Q] the ions are coupled using the motional ground state, but it was shown later that 
this requirement can be relaxed Eo|. 
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The free Hamiltonian, Hq, is given by 



N N 
n— 1 j n— 1 /e m / 

These four terms represent the energies of the excitons, of the QD phonons, of the hnear support phonons, and of the 
external electromagnetic field, respectively. Here n is the QD index, \'^j)n is an exciton eigenstate in the n'^ QD, 6„fc 
is an annihilation operator of the A;"^ phonon mode of the n^^ QD, am is an annihilation operator of the m'^ linear 
support phonon mode, and q is an annihilation operator of the l^^ mode of the quantized external electromagnetic 
field. The phonon frequencies of the support are denoted by lj*, and those of the quantum dot by u'^ . 
The coupling Hamiltonian, Hc^ is given by 

Hc^Y^ Pnokl\^j)n{'^M + E 7«r*fc|*j)n(**|^L + ^ ^nklhnkc\ + H.C. (2) 
njkl njkl njkl 

The first term is responsible for radiative decay of exciton states. The most important radiative decay pathway is the 
recombination of the electron and the hole. The second term describes the exciton-phonon interaction and gives rise 
to both pure dephasing and to non-radiative transitions between exciton states. The third term is a coupling between 
the QD phonons and the electromagnetic field. 

The interaction Hamiltonian Hj describes the coupling of the excitons to single- mode plane- wave lasers in a standing 
wave configuration. We treat the laser fields semi-classically. A QD has a permanent dipole moment due to the different 
average spatial locations of the electron and hole. In the dipole limit, we can write 

i?/ = D • E 

,(*.|e(r: - r5J)|*,)„|*,)„(vI/,|] • [euEusin(k • r?^ + (/.,)cosKt - ■ (3) 



Here r" and r)J are the position vectors of the electron and hole in the n'^ QD respectively; r"j^ is the center of 
mass location of this QD; , Ek and z^k are, respectively, the polarization, electric field amplitude, and frequency 
associated with the field mode k; (j)^ and (pt are the spatial and temporal phases of the field. The dipole limit is 
valid here since a typical energy scale for single-particle electronic excitations in QDs is 0.1 — leV, corresponding to 
wavelengths 1/fc ~ 0.1 — Ifim. For a typical dot radius i? < 5 nm, the electric field is then almost homogeneous over 
the dot. In analogy to ion trap schemes 0, the center of mass of the rt*'' QD, r^^^, is decomposed into its constituent 
phonon modes, 



•■cm ~ / ^ <-mnHm 



where are normal modes and qom is the zero-point displacement for the m}^ normal mode, qom —y/h/2MuJm, 
where M is the mass of the mode and Wm is the mode frequency. For low phonon occupation numbers, where the 
motion of the center of mass of the QD is small compared to the wavelength of the light, the Lamb-Dicke regime is 
obtained, i.e., 

Vmnk = k • C„i„qom < 1. (5) 

Therefore, we can expand Hj to first order in the Lamb-Dicke parameter 77: 

Hi = 2h'Y g^^"|*i)„(^'j|cos(z^ki - (t>t) sin^a; + ^ 77m„k(aL + am)cos0x ■ (6) 

kijn \ m / 

Here 

9^ = n{^^\^e^■ir:-rl)\^J)n (7) 

is the resulting coupling parameter between the carrier states in the n}^ QD. The second term in Eq. (^) transfers 
momentum from the laser field to the QD, thereby exciting phonon modes of the linear support. This term allows 
us to perform two-qubit operations, as described below. Manipulation of the spatial phase (p^ allows us to selectively 
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excite either the carrier transition, i.e., a change in the internal degrees of freedom of the QD without changing the 
vibrational state of the support, or a side band transition in which the internal degrees of freedom of both the QD 
and the vibrational mode of the support are changed, depending on whether our QD is located at the antinode or 
node of the laser, respectively . 

Now, let = 27r/rop be the Rabi frequency of our desired quantum operations [see Section IV |, and let T be the 
temperature. Our system must then satisfy the following set of basic requirements: 

1- Top < Tree, where t^cc is the time scale for exciton recombination. Typically t^cc = 10^"^ — lO^^sec [|60| , ^ . 

2. f2 < cjj, where cof is the first harmonic of the linear support spectrum. This requirement must be met in order 
to resolve the individual support modes. 

3. ki,T < hujf. This ensures that only the ground state phonon mode is occupied. This requirement comes from 
the Cirac-Zoller ion-trap scheme [0, where the motional ground state is used as the information bus. 

4. Dephasing and population transfer due to exciton-phonon coupling must be minimized, or preferably avoided 
altogether. 

We now discuss the details of our system in light of these requirements. 



III. QUBIT AND LINEAR SUPPORT DEFINITIONS 
A. Definition of Qubits 

In analogy to the Cirac-Zoller ion trap scheme [Q, three excitonic states will be used, denoted |0)„, |1)„, and |2)„. 
Very recent advances in ion trap methodology have allowed this requirement to be reduced to only two states [p2[ (see 



also 93 1). However, for our purposes it suffices to use the more familiar three state scheme. The states |0)„ and |1)„ 
are the qubit logic states, and |2)„ is an auxiliary state that is used when performing two-qubit operations. These 
three exciton states must possess the following properties: 

1. Dark for optical recombination: This is required so that we will have long recombination lifetimes. 

2. Dark for radiative relaxation to other exciton states: This is required to prevent leakage to other exciton states. 

3. Degenerate: The energy separation is required to be smaller than the lowest energy internal phonon, in order 
to suppress nonradiative transitions between states. We wish to make transitions between vibronic eigenstates, 
rather than creating oscillating wavepackets which would dephase as they move on different potential surfaces. 
Obtaining a large amplitude for moving between vibrational eigenstates of two surfaces depends on the exis- 
tence of two features (i) large Frank-Condon overlap Q between these eigenstates, which will be the case for 
degenerate exciton potential energy surfaces; and (ii) very narrow-bandwidth laser pulses which can selectively 



address the required states. These transitions are described in Section IV A. The degeneracy will need to be 
broken in order to perform certain operations. 

In order to choose states which satisfy the above requirements, detailed calculation of the exciton wavefunctions and 
fine structure of the quantum dots is essential. We employ here the multi-band effective mass model which has been 
employed by a number of groups for calculation of the band-edge exciton fine structure in semiconductor QDs made of 
direct band gap semiconductors [^,0 . For larger nanocrystals, possessing radii R > 20A, multi-band effective mass 
theory is generally in reasonably good agreement with experiment as far as energetics are concerned [ffsf . It has been 
used extensively for CdSe nanocrystals by Efros and co-workers |^ . While the effective mass approximation (EMA) 
has known serious limitations [S^ , and has been shown not to provide quantitative results for smaller nanocrystals 
[ [73| , it nevertheless provides a convenient, analytically tractable description, with well defined quantum numbers for 
individual states, and will allow us to perform an order of magnitude assessment of the feasibility of our scheme. 

To explain the exciton state classification resulting from the multi-band EMA, it is necessary to consider a hierarchy 
of physical effects leading to an assignment of appropriate quantum numbers. These effects are, in decreasing order 
of importance: (i) quantum confinement (dot of finite radius, typically smaller than the bulk exciton radius), (ii) 
discrete lattice structure (iii) spin-orbit coupling, (iv) non-spherical nanocrystal geometry, and facetting of surfaces, 
(v) lattice anisotropy {e.g., hexagonal lattice), (vi) exchange coupling between electron and hole spin. The electron- 
hole Coulomb interaction is neglected: detailed calculations show that this may be treated perturbatively over the 
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range of nanocrystal sizes for which the EMA is accurate . These effects lead to the following set of quantum 

numbers: rig (n^) the principle electron (hole) quantum number, Je (Jh) the total electron (hole) angular momentum, 
Le {Lh) the lowest angular momentum of the electron (hole) envelope wavefunction, Se [Sh) the electron (hole) Bloch 
total angular momentum, and the total angular momentum projection 

=mK +ms, (8) 

where mx — ±1/2, ±3/2, ±5/2, .. refers to the projection of the hole total angular momentum J^, and = 
±1/2, ±3/2, ±5/2, .. is the projection of the electron total angular momentum Jg. State multiplets are classified 
by ngieJe nhLhj^, e.g. lS'i/2l-P3/2) and states within the multiplet are labeled by F^. For a II- VI semiconductor such 
as CdSe, the Bloch states for the valence band hole states possess total angular momentum Sh = 3/2, 1/2, deriving 
from coupling of the local orbital angular momentum 1 in p-orbitals, with hole spin 1/2. The corresponding Bloch 
states for the conduction band electron states have total angular momentum Sg = 1/2, deriving from coupling of 
the local orbital angular momentum in s-orbitals, with electron spin-1/2. We consider here only states within the 
band edge multiplet, for which = Lh = 0, and Se = 1/2, Sh = 3/2. Hence the total electron and hole angular 
momenta are given by Je = 1/2, Jh — 3/2, respectively and there are a total of eight states within this multiplet. It 
follows from Eq. (||) that there is one — 2 state, two Fz — 1 states, two F^ = states, two Fz ~ —1 states, and 
one Fz = —2 state in this l<5'i/2l'S'3/2 multiplet. States within a doublet are distinguished by a superscript (L or U). 
The eigenfunctions, linear absorption spectrum and selection rules for dipole transitions from the ground state to this 
lowest lying EMA multiplet are calculated in Ref. j6^. For spherical QDs the following results were found: 

• Hexagonal crystal structure: The Fz — ±2 states constitute degenerate exciton ground states. The Fz — ±2 
states and one of the Fz = states (denoted 0^) are optically dark in the dipole approximation. 

• Cubic crystal structure: The Fz = 0'^,±1^,±2 states constitute degenerate exciton ground states, and are all 
optically dark. 

We consider here explicitly a nanocrystal made from a direct band gap material with cubic crystal structure. An 
exciton wavefunction of the 1 6*1/2 1'5'3/2 multiplet, "^p^lje, rh), can be expanded in terms of products of single-particle 
wavefunctions m (""e) and 'ipy2 tjxa-^'"'*) |H'@- -'^^ order to satisfy the requirement of optically dark qubits (recall 
condition 1 above), we construct our qubit from the Fz = —2,0^ states: 

|0) = |*o^(re,r/i)) = IV'f/2,-l/2(l"e)V'f/2, + l/2(r/0> - IV'f/2, + l/2(re)V'3/2,-l/2(rft)) 

|1) = |*-2(re,rft)) = |^f/2,_l/2(l-e)V'f/2.-3/2(l-/.))- (9) 

The auxiliary level for cycling transitions is taken to be the Fz ~ 2 state: 

|2) = |*+2(re,r,0) - l^f/2, + l/2(re)V'f/2,+3/2(l-/0>- (10) 



Explicit expressions for the electron and hole wavefunctions are given in Appendix IX. The states |0), |1), |2) are 
degenerate and have equal parity [determined by (—1)^-]. 

Naturally, for nonspherical, noncubic, and/or indirect band gap materials, other states may be more appropriate. 
It is only essential that they satisfy the requirements above. In this paper we shall use primarily the EMA states 
described above for cubic nanocrystals of direct gap materials, because they illuminate in an intuitive and quantifiable 
manner the difficulties associated with our proposal. However, in the discussion of feasibility (Section ^), we will also 
present results obtained with qubit states obtained from tight binding calculations for nanocrystals constructed from 
an indirect band gap material (silicon). 



B. Quantum Linear Support 

In order to determine whether quantum computation is possible on such a system we need to examine also the 
properties of the linear support. The support is made out of K small units, e.g., unit cells or monomers. We can 
write the displacement of each unit as a sum of normal modes, 



Xfc = ^ Crnk<im = ^ Cm^qom (a|'„ + «„). (11) 
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The zero point displacements for a homogeneous support are 



qom = \/'hl {2Xluj„i), (12) 

where A is the hnear mass density and I is the length of the unit. The lowest energy modes will be long- wavelength 
transverse modes. Since the wavelengths of the modes of interest are large compared to the separation between 
neighboring units, we can approximate the support as being continuous. 

In many cases, a sparse number of attached QDs will have only a small effect on the normal modes of the support. 
The validity of this assumption depends on the materials chosen and will be discussed more thoroughly below. For 
now we will calculate all of the relevant properties assuming point-like, massless quantum dots, consistent with 
our assumption that the spacing between the dots is larger relative to their intrinsic size. For the n*'' point-like QD 
attached to unit cell k, with one dot per unit cell, we can identify the dot and cell normal modes expansion coefficients. 
We then have Cmn = Cmk, where Cmk and c„m are, respectively, the coefficients relating the displacement of the k*"^ 
unit cell and n*"^ QD to the displacement of the m}^ normal mode [Eqs. (jj) and (^l]), respectively]. For a continuous 
support, the set of Cmk becomes a function Cm{x) that is the normalized solution to the wave equation on the support. 
Any specific value of Cmn can now be written c„i(x„), where Xn is the position of the n^^ QD. 

The two most common types of linear continuous systems are strings and rods. 



1. Strings 

In a string, the resistance to transverse motion comes from an applied tension, ??. The dispersion relation for the 
frequency of a string in mode m is 

- ^^^k 

where A is the linear mass density and km is the wavenumber. The normalized solution to the transverse wave equation 
with fixed ends is given by 

where km = rrni/L, I is the unit length, and L is the string length. 



2. Rods 



In a rod, the resistance to transverse motion results from internal forces. This leads to a different dispersion relation 
and, consequently, to a different solution Cm(x). The transverse modes of a rod can be defined in terms of the length 
i, density p, the Young's modulus Y , the cross-sectional area A, and the second moment of A (or the massless moment 
of inertia of a slice of the rod), I. As shown by Nishiguchi et al. [Q, the long wavelength phonon modes (A > 1000 
A) are well described by the classical dispersion equation 

^ra - km\j 

where is the moment in the direction of the displacement. 

The transverse normal modes for a clamped rod are well known | |9^ , resulting in the solution 

Cm{x) = Nm [sin(fc„iL) - sinh(fc,„i)] [cos(fcmx) - cosh{kmx)] 
— [cos{kmx) — cosh(fcma;)] [sin(fc,„L) — sinh(fc„jL)] . 

Here Nm is a normalization constant proportional to ^J\JL. The values of fc„i are not known analytically, but can 
be shown to be proportional to XjL. 
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C. Approximations 



The important parameters characterizing the support are ujI, the frequency of its first harmonic, and the product 



This product is the quantum dot displacement resulting from the zero point motion of mode m of the support. We 
shall refer to it as the dot modal displacement. The above discussion of vibrations in the support assumed massless 
QDs, motivated by the assumption that they have negligible spatial extent relative to the distance between them. To 
investigate the effect of the finite mass of the quantum dots, we computed numerical solutions of the coupled vibration 
equations for strings and rods having finite mass increments located at discrete points, simulating the attachment 
of finite mass quantum dots. These numerical calculations show that for sparsely spaced dots of mass small enough 
that the total weight is the same order of magnitude as the weight of the support alone, the resulting value of Snm 
remains unaffected to within a factor of 2 by the presence of the dots (see Figure A simple way to approximate 
the presence of the QDs and retain an analytic solution is then to replace the linear density of the support by the 
average combined linear density of QD and support. 

Since we are interested here in order of magnitude estimates of feasibility, we will approximate Cm(xn) by y^2l/L. 
This approximation corresponds to the maximum value of Cm{xn) for a string, and to approximately the maximum 
value of Cm(x„) for a rod. Since the larger the dot displacement, the larger the coupling between dots, this means 
that our estimations of number of operations will be an upper bound. 

These approximations combined with Eq. yield the following equation for the Lamb-Dicke parameter. 



r]m'nk2 = k2 ■ Snm = k2Snm COS 9 = k2\ — COS 6*, (14) 

V -IVlUJm 

where 9 is the angle between the modal displacement and the direction of the laser beam, and M is the total support 



mass: M = LA. Note the inverse power dependence on M in Eq. (|14D. We shall see in Section VB that the massive 
nature of the linear support and the resulting small value of the Lamb-Dicke parameter provides the major limitation 
for our system. 



IV. QUBIT OPERATIONS 
A. One-qubit Operations: Coupling of Dots to Light 

1. Derivation of the Interaction Hamiltonian in the Rotating Wave Approximation 

While dipole transitions in QDs are similar in principle to dipole transitions in atomic systems, the strong coupling 
to internal phonon modes adds an additional complexity. Consider the modifications of Eqs. (|l|)-(^) for a single QD 
unattached to a linear support, interacting with a single laser field, with the QD located at the anti-node of the field, 
i.e., with the sin(k • r'^^ + (p^) term in Eq. (||) vanishing. Thus, omitting the linear support term: 

j k ijk ij 

We separate this Hamiltonian into two parts, and Hi : 

Ho^Y. KI*.^-)(*.^-| + E ^""'bibk + E Hjfcl*.7->(*.l(4 + bk) (15) 

j k jk 

Hi - E H-^fel*.>(*'^l(&I + bk)+Yl 2%k'l*^>(*.l cos(i^kt - </>*) (16) 

i^j.k ij 

The "free" Hamiltonian Hq may be diagonalized by a displacement transformation. Let Dk{a) be the unitary dis- 
placement operator: 

Dk{a) EE e""!-"*^- = Dki-ay. 
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The displaced phonoii operator djk is defined as 

djk = Dk{~ajk)hkDk{ajk) = bk + ajk, 
and satisfies standard boson commutation relations: 

[djk,dl,k,] = Skk' 
[djk,dj>k'] = [4fe'4''c'] ^ 0- 

Note that for real ajk we have ajk{b\ + bk) — dji^djk — ctjk^ — b\bk- Letting Ujk = ^jjkl<^k a-nd inserting a complete 
set of exciton states into Eq. (p^, we find: 

H^^Y. Ki*^) + E E 1*^) (^I^fe - '^^^^^l + ^^)) 



The eigenstates of d^j^d^k are labeled \njk) , where c?]fedjfc|njfe) = n\njki and {njk\mj,k') = hk'F^mk^ The i^^^^, are 
Franck-Condon factors M], describing the overlap of vibrational eigenstates between different excitonic states j and 



j' . We can then rewrite Hq as 



i/o = Y.n^'j+\^,){^,\ + Y.hu;t\^,){^,\dy,k, 

j jk 



where uj'j^ ^ ujj — (jjjk/^k)^ is the renormalized electronic energy level. 

We transform to the interaction pictm'e defined by Hq: Hj = exp {iHot/Ti) Hj exp {—iHot/h). To do so, it is useful 
to insert into this expression two complete sets of displaced oscillator states belonging to different excitonic states i 
and j: 

h X Ij = 0E l"'^)^"*^! ^ 0E = 0E l"''=^*^'"j'=l-^"™fc- 

k n k' m k nm 

Changing variables from bk to djk in Eq. ([l^), and transforming to the interaction picture now yields, after some 
standard algebra: 

ffj — ^iHot/hjj^^-iHot/h 

= E ?^7,.fe|*.>(*,|e'<*(g)E^™^l"^fc)K-fc|e'"''""'"'* 

i^j,l k^l n,m 

111 \ 

+E^3ki*.>(*.k-"^^''''*~"^^*'^'^(8)E^-^^ (17) 

where = 0;,''+ - 

While this expression appears very complicated, it can be drastically simplified under certain reasonable assump- 
tions. First, note that for single-qubit operations we need to consider only two exciton states j'l'a) and l^*;,). The 
first term in Hj essentially describes non-radiative transitions between exciton states due to phonon emission. Under 
the assumption that the phonon modes are initially unoccupied, we can choose the states j'l'a) and \^/f,) such that 
they have a negligible propensity for nonradiative transitions, i.e., they are protected against single phonon emission 



(recall condition 3. for "good" qubits in Section illA). This means that we can effectively set all 7jife to zero, thus 



Note that n refers here to the occupation quantum number of the internal phonon modes, not to the quantum dot index. 
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eliminating the first term in Hj. This important simplification is treated in detail in Section VA below. Thus we are 
left with: 



Hi = ngf |vI/,)(vI/,|e^**-'*('^>'— (g)^F„t,|n„fe)(m6fc|e'"^("-™)* + H.c. 

k n,m 

We then tune our laser on resonance such that Vk — ^at^ ^^'^ make sure that the laser spectral width is much smaller 
than the lowest quantum dot phonon frequency, ujf. This allows us to make the rotating wave approximation (RWA), 
i.e., eliminate all terms which rotate faster than ojf, which leads to 

Hi = hg^'\^a){^b\e''^*<S)T.Pnnk\nak){nbk\+il.c. (18) 

k n 

This RWA interaction Hamiltonian, Eq. (^8|), is very similar to the familiar two- level system Hamiltonian used 
extensively in atomic optics |p8|] . However, the strength of the interaction is modulated here by the Franck-Condon 
factors, F!^^f.. To allow the simplification of the Hamiltonian from Eq. ([l7| ) to Eq. ( |l8| ) requires a judicious choices of 
laser intensities and states. In our scheme the occupation n of all phonon modes will be initially zero. Using Eq. (^, 
it is useful to then introduce the factor 

= gl'XlFSL - \/^^n^o°ofe(*a|ek • (re - r,)|*,), (19) 

k k 

which corresponds to the Rabi frequency for an on-resonant transition. Ik is the laser intensity and a — e^/ (47r£o^c) 
is the fine structure constant. 



2. Raman Transitions 



Since we wish to use near degenerate states of equal parity for our qubits, we cannot employ dipole transitions. 
Hence we use Raman transitions. These connect states of equal parity via a virtual transition to a state with opposite 
parity. Recall that parity is determined by (—1)-'^^. Suppose we start in the |1) = |V'f/2 -i/2^^<^)'^i/2 state, 
for which — —2. We can then make transitions through a virtual level \v) that has opposite parity (e.g., Fz = ±1), 
to the state |0) having F^ = 0. Figure |^ provides a schematic of the coupled QD-laser field system, showing the 
levels |0), |1) and \v) together with the fields required to cause a Raman transition. Under the assumptions that only 
two laser field modes ki and k2 are applied, and in the rotating-wave approximation, the standard theory of Raman 
transitions | ]99| leads to the following expression for the Raman- Rabi frequency between an initial state \i) and a final 
state I/) : 

^'^ti - (20) 

Here ki and k2 are chosen such that the detuning, IS. — Uj ~ Ui — i^ci — ~'^f ~ ^2 1 J the index of an intermediate 
exciton state chosen to provide a minimum value of A. For single qubit transitions, both lasers are aligned such that 
the QD is positioned at antinodes. 

For QDs possessing cubic crystal structure and composed of direct band gap materials, we have found it advan- 
tageous to use the F^ — 1 and F^ = —1 states of the 15*1/21^5/2 multiplet as the intermediate state. An exciton 
wavefunction of the l>S'i/2lP5/2 multiplet ^E*^ (feT'i'h) can be expanded in terms of products of single-particle wave- 
functions "03/2 m (^e) ^^'^ "^5/2 niK^'''^^' P^ i Poj The intermediate F^ = ±1 virtual states can be written: 

IV'f/2,=Fl/2(l"e)V'f/2,±3/2(r/i)) + |l/'f/2,±i/2 (re ) V'f/2,± 1/2 (l"/i)> 

The Raman-Rabi Frequency, riRaman, can then be adjusted by increasing the electric field intensity and by reducing 
the detuning from the intermediate level. We will describe in detail in Section |VB| what range of values of intensity 
and detuning are allowed. 
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B. Two-qubit Operations: Coupling Quantum Dots, Quantum Supports and Light 



Our two-qubit operations are equivalent to those of the Cirac-Zoller scheme . The use of optical Raman transitions 
to implement this scheme has been extensively explored. |^ In our case, we apply the Hamiltonian of Eq. with 
two lasers ki and /c2, of frequency vi and 1^2 respectively. For two-qubit operations, the quantum dot is centered at 
an antinode of fci and at a node of k2. Switching to the interaction picture and calculating second-order transition 
probabilities to first order in r/, one obtains the following effective Hamiltonian: 



H. 



nfi 



H.c. 



(21) 



Note that the nodal and antinodal lasers result in an effective Hamiltonian in which 77 depends only on the nodal 
laser ^2 • This differs from the effective Hamiltonian derived for Raman transitions when travelling waves are used ||^ . 
The lasers are chosen to have a net red detuning, uj f — tOi — vi + V2 — —ujm'- In the RWA (i.e., eliminating all terms 
rotating at 2ujm), with (j)2 — (pi = this yields 



h: 



nfi 
eff 



/ik2ki 



H.c. 



(22) 



where 



nm' 1^ 



Q/ik2ki 
2 ' 'Raman • 



(23) 



This combined QD-linear support operation transfers the n^^ QD from state i to /, with an accompanying change 
of one quantum in the phonon mode m' of the support. A schematic representation of this operation for the qubit 
states \ipi) = |0) and IV'/) — |1) is shown in Figure ^. Choosing interaction times such that t = fc7r/(2r2^^?''^) where 

k is an integer specifying the pulse duration, we can write the unitary operator exp (^—j-H"^h^ — U^^^it) as 



nfi 



exp 



(24) 



In order for the Cirac-Zoller scheme to be successful, the phonon mode of interest, to, must start with zero occupation. 
The applied operations take advantage of the fact that the zero occupation phonon state is annihilated by the lowering 
operator: a|0) = 0. The sequence of unitary operations J7c-phaso = C/f ^°[/^ 20jjnio f^^^Qj^ results in a controUed-phase 
operation between quantum dots n and n' , i.e., it causes the second qubit n' to gain a phase of —1 if the first qubit 
n is in the |1) state, and no additional phase if the first qubit is in |0). This is equivalent to the matrix operator 
/ — 2|1)„/|1)„ „/(l|„(l|. The time required to perform t/c-phase is then 2 x for ion n plus 1 x tt for ion n' , i.e.. 



_ — /ctnn'\-l _ //-,/ik2ki N -1 . /(-,/ik2kiN-l 



(25) 



Since the ions n and n' are identical we will use the approximation that ^^^f^^ ~ ^n'm'''^ the remainder of this 
work, hence 



2tt "™ 2n ^ 



(26) 



The inverse of the average rate " can then be taken as a measure of the gate time, i.e. 
controlled phase operation. We define ^l2 as the sideband interaction strength, 



o o/^k 
"2 — f7rTmk2"Rai] 



ki 

Raman " 



of the time for the two-qubit 



(27) 



Calculation of ^{^^^^^^ was described above in Section IV A [Eq. (p^)]. We can obtain the Lamb-Dicke parameter 
Vm'nk.-2 from the decomposition in Eq. (p^. This now allows specific evaluation of the contribution from the linear 
support to the Lamb-Dicke parameter 77. As described in Section III C , we approximate 77 as being independent of 
the specific dot and from now on will drop the dot index n. 
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C. Input and Output 



Since the qubit states do not include the ground state of the quantum dot, initiahzation wiU generaUy require a 
transformation from the ground state of no exciton to the defined qubit state. This can be accompUshed by applying 
magnetic fields which will mix dark and light states allowing for optical transitions. If the magnetic field is then 
adiabatically removed, one is left with population in the dark exciton state only. Qubit measurements can be made 
by using a cycling transition, in analogy to ion traps . 

We conclude this section by summarizing in Fig. the relative energy scales involved in our proposal. 



V. FEASIBILITY OF QUANTUM LOGIC 

In this section we address in detail the question of the limitations imposed on our system by various physical 
constraints. We start by considering the issue of decoherence due to coupling of excitons to the internal nanocrystal 
phonon modes, and propose a solution to this problem. We then study the issues of scaling arising from the trade-off 
between massiveness of the support, laser intensity, and the need to maintain a large ratio of operations to exciton 
recombination time. We find that the allowed size of our proposed quantum computer depends on the assumed 
threshold for fault-tolerant computation. 



A. Decoherence 



According to current analysis of experiments on nanocrystal quantum dots |100|, exciton dephasing derives pre- 
dominantly from the diagonal phonon exciton coupling term of Eq. (H): 

J2\^)(^\("f*^A+"fM- (28) 

Here jkj is the self-coupling of an exciton state j via phonon k, and bk is the lowering operator for the A;*'' phonon 
mode in the ground electronic state. In the ground state the coupling is zero, and all excited electronic states have 
potential energy surfaces which are shifted with respect to this ground state. We desire to eliminate dephasing due 
to the first-order phonon exciton interaction. In the typical experimental situations in which dephasing has been 



studied in the past, dephasing occurs on a timescale of nanoseconds for small dots |10C]. This rate is extremely rapid 
compared to the experimental recombination time of dark states 10~^ s for direct band gap materials such as 
CdSe [Q). The reason for such fast dephasing is twofold. First, the vibrational stationary state of the first electronic 
level becomes a moving vibrational wave packet on the upper electronic surface, because the spectral width of the 
pulse is too broad to distinguish vibrational eigenstates. Second, the QD is embedded in a solid state medium where 
the vibrations of the nanocrystal are then coupled to vibrations of the larger lattice. The phonons of the QD can 
be treated as analogous to damped cavity modes in atomic optics [Q. In the case of strong coupling, one finds 
from numerical simulation that the dephasing between any two states j and j' is related to the rate of phonon mode 



excitation. The latter is proportional to — Tfcj'P for each mode k [100|. This conclusion of f ast d ephasing agrees 



with the analogous result for a leaky optical cavity |£8[ as well as with the results of experimental 1 101 ] and theoretical 



I IOC I analysis for embedded semiconductor nanocrystals. 

The dephasing can be reduced in three ways. First, the coupling of the QD phonon modes to external phonon 
and photon modes can be reduced by judicious choice of nanocrystal geometry and material. In our case the QD 
can dissipate phonon modes only to the support. In the limit of no coupling to external modes, there will be no 
dephasing but the time required for recurrences could limit our gate repetition rate. Although the oscillations will be 
fast, the oscillations for different phonon energies will be incommensurate with one another. This could introduce a 
slow quantum beating between ground and excited electronic states, which would have the undesirable consequence 
of requiring gate durations to equal a full beat cycle. 

Second, one can find a set of electronic states \ j) such that 

Ikj - Ik]' = Vj,/, and Vfc. (29) 
Physically this condition represents a set of electronic states which create the same potential energy surface for nuclear 



motion. This elimination of decoherence by degeneracy is an example of a decoherence-free subspace 1 102 -104]. The 
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Jahn- Teller effect implies that no two such states should exist, because there will always be a phonon mode which 
will distinguish between these states due to non-linearity |105|. However, in the linear approximation we have 



Ikj = (j|7fe(r,i) +7fe(re)|j). (30) 
The deformation potential coupling operator, 7fe(r), is a function of the phonon modes and is expressed as 

7fc(i") = 7nim(r) = -EdV • u„;„(r), (31) 
where Ed is the deformation potential, and Unim{f) is the coordinate representation of the normalized spheroidal 



phonon mode of level n with angular momentum I and projection ni. Following Takagahara [100|, the spheroidal 
modes can be written as 



U„im(r) =W- {pnlLlra[hnir) + Qnl'N lrn{knir)) , (32) 

where p is the nanocrystal density, ujnim is the frequency of the spherical phonon nlm, Lijn{hr) — ■iV^'/m(ft-r), 
Nim{kr) — iy X V X r'^im{kr), and '^{kr) = ji{kr)Yj^ {^) . ji (r) is an order spherical Bessel function, Y/"-{^) is 
a spherical harmonic, fc„ and /i„ satisfy stress- free boundary conditions at the surface, and Pni, Qni are determined by 
normalization. One can then write 



7«im(r) = -Ed\ —-^ — Pnihniji{hnir)Yi™-{Vl). (33) 

For a cubic, direct gap nanocrystal, the states \j) are states of well defined angular momentum projection. Since 
the y^™ in Eq. ( p}3[) connects states with equal projection, the only phonon modes which can have non-zero matrix 



elements in Eq. (pO[) are those with m ~ [106|. The resulting matrix elements will be independent of the sign of 
the exciton angular momentum projection Fz [Eq. (|^)], i.e., the to = phonon modes cannot distinguish between 
exciton states having or —F^. Therefore, in the linear approximation, the states|5'_2(^ej ^h)) and \'^2ifeii"h)) will 
not dephase with respect to each other. Recall that we took these states as our qubit state |1) and auxilary state |2) 
states [Eqs. (|),(|l^)]. 

Third, and most importantly, one can change the way in which transitions are made. In the above two situations, the 
motional wave packet of one electronic surface is transferred to another electronic surface without changing shape, i.e., 
the Franck-Condon approximation holds. However, such a transition requires either a broad laser or a fast excitation. 
This is not actually the regime of relevance here. The scheme outlined in this work requires selective excitations of 
sidebands whose energy separation is orders of magnitude smaller than the quantum dot phonon energies (see Fig. H). 
Therefore, we will be performing transitions from one vibrational eigenstate to another vibrational eigenstate. Such 



transitions were described in Section [V with respect to the ground vibrational state. Consequently, the scheme 
proposed in this work is not affected by fast phonon dephasing. 

Exciton states recombine and thus decay to the ground state by both spontaneous emission of photons and phonons. 
We denote the recombination lifetime r^e. Nanocrystals are known to have dark state recombination times ranging 
from 10~^s to 10~'^s, depending on the material chosen ||60|,|73|,|9l|,p5[ . In our system, Tj-e will be the fastest decoherence 



time for individual qubits. One could potentially suppress radiative recombination by placing the whole system in 



a cavity |107]. Classical calculations of Roukes and co-workers show that nanoscale rods at low temperature have 
high Q values: Q > 10^" |108|. This implies that the rods are only very weakly coupled to their environments and we 



can therefore assume that in the quantum regime, the high Q will lead to favorably long decoherence times. Another 
possible source of decoherence is laser scattering from the support (as opposed to the QDs). The magnitude of 
scattering is dependent on the difference between the spectra of the quantum dots and the electronic and vibrational 
modes of the support. A detailed analysis of this potentially important decoherence mechanism is beyond the scope 
of this paper, due to the many possible materials available for both linear supports and quantum dots. Ideally, one 
would like to choose a support which has an optical window at the frequencies of the lasers used to perform the qubit 
operations. Estimations based on the Raman transitions proposed above suggest that this optical window needs to 
be between 0.1-1 eV and possess a minimal width of 0.1 meV. Conversely, one can take advantage of the optical 
tunability of the QDs via their size to construct a QD of such size that its transitions are compatible with a specific 
optical window suggested by the material properties of the support. 
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B. Parameter Space 



We now explore for what range of physical parameters quantum computation is possible within our proposed 
scheme, by estimating the two-qubit gate fidelity, T . This fidelity is defined as the trace overlap between the desired 
final and the achieved final state: T — minpf, TrA/9o^^5(po)j where A is the exact unitary operator for the gate, B 
is a superoperator describing the actual evolution of the system which takes the initial density matrix po to a final 



density matrix p/, and po ranges over all possible input states. For our two-qubit operations described in Section |1VB 
the fundamental operation is the population transfer to the red side band, A — J7"^'^[Eq. (p^]. B describes both 
the unitary evolution caused by application of the lasers and the decoherence due to loss of quantum information 
to the environment. Note that even without decoherence and unknown laser noise, T can still be less than unity, 
due to deviations from the approximations used to derive A. Most importantly, deviations from the rotating wave 
approximation can lead to unwanted population in spectator states. 

The resulting value of T is determined by two constraints: the decoherence time of the system and the spectral 
resolution of the gate. As a result of the use of the phonon bus in the two qubit gate construction, both this proposal 
and the ion trap proposals have gate times that are dependent on the number of qubits. Therefore the quantity 

of interest is the fidelity of a sideband operation on an array of iV quantum dots. First, we note that J- is limited by 
the recombination time of the qubit states, Tre = l/F. Naturally, the recombination time Tic must be larger than the 
sideband operation time, ta — 7r/(2r22), if the operation is to be successful. One can then define an upper limit on 
the fidelity that takes into account the statistically independent recombination of the exciton states of all N quantum 
dots. This background fidelity assumes that Ti-e is the same for all N quantum dots and it does not account for errors 
deriving from the interactions with the driving laser field. The background fidelity is then 

T is also limited by the spectral resolution. We assume that the difference frequency between the lasers is tuned to 
be resonant with the energy difference between the two states of interest. Figure ^. The energy difference between the 
states includes the relative Stokes shift induced by both the lasers and the internal phonons. Omitting the adjustment 
of the laser frequency for the Stokes shifts will lead to gates of reduced fidelity | 109| , 110 |. The primary concern is that 



the phonon modes be spectrally resolvable. As described in Section ||, the use of nodal and antinodal lasers allows 
us to transfer population to the vibrational sideband and at the same time forbid population transfer to the carrier, 
i.e., to the excitonic states of the QD. This constrains the operation frequency ^2 to be smaller than the separation 
between phonon modes, Acj = — uj!^^. For small m, and for Ato « ajf, we can then write this constraint as 

n2 < col (35) 

The population transfer to the off- resonant state will be of the order of ^gi_^^gyi , the result for a two level system 
interacting with a periodic perturbation of strength g that is off-resonant by a frequency difference 6 |9^]. In our 



case, g = ^2 and 5 — Alo « wf. However, numerical calculations [ 109 | have shown that the population transfer 

2 r 

to the off- resonant state for an ion trap system is more accurately estimated by 2 ^^^^^.^^^ . Note that in Ref. [ |109 
g=r2/2, the coupling to the carrier transition, and 6 = Uz, the ion trap mode vibrational frequency. The advantage 
of the standing wave laser configuration is now apparent. For the travelling wave laser configuration, the population 

transfer to the off-resonant state is 4(^^^^j^yi!l(^^syi ~ 2 (^jj^^ (or, in the notation of Ref. [109|, ^ (l7") )• However, 

in the standing wave configuration, one finds that the population transfer is '2j(ji^^^^f(jjTyi ~ 2 . Since 77 <C 1 

the off-resonant population transfer is significantly reduced when one uses the standing wave configuration. We can 
now write down a fidelity which takes into account both the background fidelity, Eq. (Q), and the population loss to 
the most significant spectator state. The fidelity per sideband operation A is then 

ttNT ./1^2^^ 



in the standing wave configuration. We emphasize again that if the laser fields are used in a traveling wave configura- 
tion, the fidelity is significantly decreased due to transitions to the carrier state, resulting in ~ 1 — — 

One can now maximize the fidelity with respect to the coupling strength fl2 for a sideband operation made on an 
iV qubit array. Since 
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nr-=(^^)"\ (37) 

for given T and Lof, the maximum fidelity can be written as 

Evaluation of the optimal operation frequency depends then only on the factors in Eq. (^). Inspection of the 
contributions to ^2 [Eq. (p7|)] shows that the underlying adjustable parameters controlling the fidelity in general, 
Eq. (p6|), are the intensities /i, I2, the detuning A, and the Lamb-Dicke parameter 77. However, the schematic shown 
in Figure |4| shows that there are some additional constraints. Thus, it is essential that the inequality cul < A < ujf 
be satisfied in order to avoid unwanted coupling to both the internal phonons and the linear support phonons. In 
addition, we require that both jfi^*^ | and 77oik2 l^Ll I ^^^^ smaller than A, in order to avoid populating the intermediate 
level, I^Pj). Hence, the internal phonon energies define an energy scale which also constrains our system (see Fig. 
Furthermore, combining Eqs. (pO|), (p7|), and (|37|), one finds that 



\K^oi^^< I - ^ 8 j • (39) 

Analysis of three level systems has shown that in order to maximize Raman population transfer between states 
and l^*/), the coupling strength between j^fj) and j^'^) and l"^ j) and j^Pj) should be equal M]. In the atomic case this 
usually implies that the respective Rabi frequencies between electronic states, |il^'J and |^^k2li ^-^e equal. However, in 
our case, with the use of nodal and anti-nodal lasers and coupling to the support phonons, the equivalent condition 
is that 

m^voi^j^i^j- (40) 

Therefore, manipulation of Eqs. (^9|), (|39| ) and (^o|), allows us to determine the intensity values, /i and I2, necessary 
for maximum fidelity operations: 

HAiM^^Y' ^ (41) 

hA /7r(wf)2iVry^^ 1 1 



/2 ^ _ (42) 

2™V 8 ; '^oik. |(*/|62-(re-r.)|vI/,)n,^^/o;P 
_ A\L firiujlfNTY^^ 1 1 



2na\ 8 J \^2\' \{^ f\e2 ■ {ve ~ r,)\^ ,) Ui F^^il' 



Equations (38) , (|4l|), and ( |42| ) summarize the limits to implementation of this quantum dot-quantum linear support 
scheme. To maximize the fidelity we need to increase the frequency of the phonon bus, Luf. However, as this 
frequency increases, the increased intensity necessary to reach the maximum fidelity will lead to unwanted evolutions 
not considered in our simple fidelity equation, Eq. (^8|). These unwanted evolutions include quadrupolar excitation 
to higher electronic states. Such transitions will not be removed by the use of nodal and anti- nodal lasers 
Consequently, it is useful to define maximal laser intensities, I™^^ and I^^^, such that Eq.(|38|) is valid for Ii < I^'^^ 
and I2 < I™^^. Therefore the magnitude of is restricted and thereby imposes a constraint on the phonon spectrum 
of the linear support. Furthermore, scalability is also limited by the additional unwanted evolutions, since the required 
intensity to achieve maximum fidelity also increases with N . 

There are also physical constraints on the density and length of the support. We assume that the minimal linear 
density would be provided by a chain of carbon atoms, for which we estimate Ao = 10 amu /A. The length of the 
support, L, is determined by the number of QDs, N , and by the spatial width of the laser, I. Thus for identical dots, 
we have L = IN. One could use QDs having no spectral overlap, obtained by making the dots of sufficiently different 
sizes, in order to achieve more qubits per unit length. 

Notice that when one rewrites L in terms of N that the intensity of the nodal laser, I2, has a stronger N dependence 
than the intensity of the anti-nodal laser, Ii. Physically, this is due to the increased inertia of the system and is 



quantified by the Lamb-Dicke parameter rj. From Eq. (H) and Section (HIE), 77 cx M where M is the total mass 
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of the support, M = XL, so that I2 ck. M [Eq.(^2|)]. This coupled with our expression for the maximal gate fidelity, 
Eq.( |38|), y ields I2 oc N^/^ . Note that we have used the Lamb-Dicke parameter consistent with the definition made in 



Ref. 1 109 1 which has an inverse N dependence. 

To determine the scalability of our system, we examine the maximum number of QDs which can be sustained by a 
support having given values of and A/Amin, and provide an acceptable level of fidelity for the sideband operation, 
A ~ [/"^°. We do this by requiring the following three constraints to be simultaneously satisfied: i) T = J-max, ii) 
I2 < I 2^^^, and iii) jTmax > 1 — £, where e can be thought of as the error rate per gate frequency. The first condition 
states that maximum gate fidelity, Eq. (^8|), is achieved given the support and quantum dot parameters iof and F, 
respectively. The second condition states that the laser frequency I2 should not exceed the maximum allowed value 
(see above). Eq. (|4^) together with the considerations in the previous paragraph shows that I2 is dependent on the 
number of quantum dots N. For the range of parameters considered here, Ii is always smaller that I^'^^. Hence the 
maximum number of qubits will be determined by the intensity threshold of the system at a node of the laser field. 
The third condition ensures that one is able to perform an operations with success greater than a certain threshold 
value (equal to 1 — e). Combining these inequalities leads to limits on the number of qubits for a given system. 
Conditions i) and iii) can be manipulated to yield the following constraint on N: 

On the other hand, conditions i) and ii) yield a constraint with an inverse power dependence on tof. One finds that: 

1/4 / \ -^/^ 

N < {u;tr'/\iri'^' (^^^ |^|^|k,n(vl.,|e,.(re-r,)|vI/,)n^o'JzPj (44) 

One can then analyze iVmaxj the maximum allowed value of N, as a function of the linear support frequency . The 
combination of Eqs. ( ^ ) and ( ^ ) results in a cusped function for iVmax and is discussed in detail in Section VI for 
both a direct band gap semiconductor (CdTe) and an indirect band gap semiconductor (Si). 

The above discussion has focused on the fidelity for a single component operation, A — U^^^, of the C-phase gate 
f/c-phaso = Ui^^U2 ^"J/f^*^. We have termed this a sideband operation fidelity. Similar arguments may be made to 
derive the full C-phase gate fidelity, resulting in the expression 



This is lower than the sideband operation fidelity, both because of the effect of multiple couplings to spectator states 
and because of an increased operation duration (rc-phase = 4rA). In the presentation of numerical results in the next 
section we shall refer only to the prototypical sideband operation fidelity, T of Eq. (p6h . 



VI. NUMERICAL ESTIMATES FOR SPECIFIC NANOCRYSTAL SYSTEMS 

A. CdTe 

CdTe nanocrystals are an example of direct band gap cubic crystal semiconductors QDs. Using parameters found 



in Landolt-Bornstein | 111 |, we have p erfo rmed the calculations summarized in the previous sections, using previous 



results of Efros ||60[ and Takagahara [100| for the EMA analysis. Although the dark states have infinite lifetimes in 



the EMA approximation, both experimental [112| and tight binding calculations for the analogous CdSe system [ [73[ 
yield radiative recombination rates F ~ 10^. 

We analyze here nanocrystals with i? = 20 A. For this size, the frequency of the lowest internal phonon is ujf = 
2.45 X lO^^Hz. Assuming a minimal separation of a single order of magnitude between the energy spacings uj'^ 
and A, we set A = lO^-'^Hz. The energy separation between the IS1/2IS3/2 and IS1/2IP5/2 multiplets is ~ 0.4 
eV in the EMA, which leads to required wave vectors fci ~ k2 ~ 2.1/i m"-'^ for the irradiating lasers in the two- 



qubit gates. For the specific CdTe states introduced above (Section III A ), we calculate the dipole moments to be 
(*!_Y(i-e,r/,)|e2 • r|0) = 0.11 R , (re, rft)|ei • r|l) = -0.013 R where £2 = ^{x + ty) and ei = ^(x ~ ly). 

Furthermore we have calculated the Frank Condon overlap to be Yli Foot'^ = 0.98 and Yli Fom^^ — 0.98. We assume 
that the spatial width of our lasers is diffraction limited. A reasonable estimate of this width is then I = 3 /xm. 
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At constant frequency, an increase in the number of qubits requires an increase in the laser intensity in order to 
maintain maximum fidehty operations. We estimate that I^^^^ = lO^^W/cm ^ is the intensity at which the nonresonant 
quadrupole interactions begin to rise in CdTe quantum dots. However, the intensity could potentially have stricter 
limitations depending on the spectra of the specific support chosen. As mentioned above for the range of parameters 
we have examined, the intensity of the anti-nodal laser is weak enough that it does not lead to unwanted time 
evolutions. 

In Figure |[ iVmax is plotted as a function of ajf for two linear support densities and for a modest threshold of one 
error every ten operations (e = 0.1). At low frequencies, Eq ( ^ ) limits A'^max and increasing the values of ujf leads to 
larger values of A^ for a fixed J-max- In contrast, higher frequencies require stronger laser intensities [Eq.(^] so that 
eventually the limits on the intensity given in condition ii) begin to reduce the maximum possible number of quantum 
dots, leading to the turnover in Figure ^. Figure Halso shows that the optimal value of A^max, which we denote by Nc 
, is reduced for larger support densities. In Figure^, we now plot Nc as a function of the error threshold e, for a range 
of linear densities A. We see that even for two qubit quantum devices one must allow e > 0.02, or approximately 
one error every 50 operations. Even at the modest threshold value, e = 0.1, one can only support 7 qubits. Clearly, 
CdTe excitons are thus not a good candidate for scalable qubits within this scheme. The underlying reason is that 
the recombination time of the dark states, while longer than the operation time, is not sufficiently long to provide 
high fidelity operations. 



B. Si 



Si and other indirect band gap bulk materials exhibit longer exciton recombination lifetimes than direct band gap 
materials such as CdTe. Although EMA descriptions of Si nanocrystals exist, many subtleties are required to obtain 
accurate excitonic states These have also been calculated in semi-empirical tight binding approaches as 

well as via pseudopotential methods ||7^. The advantage of tight-binding descriptions is that the optical properties 
of the nanocrystal can be determined with inclusion of realistic surface effects [113|. We estimate the feasibility of 
using Si nanocrystals here using the detailed Si excitonic band structure previously calculated within a semi-empirical 
description | ]95[ | . In order to suppress phonon emission we choose states which correspond to cither the exciton ground 
state, or lie within the minimal phonon energy of the exciton ground state. The minimal phonon energies ujq are 
taken from EMA calculations made by Takagahara |10C], and are approximately equal to 5 meV for a nanocrystal 
of 20 A radius. One disadvantage of the tight-binding description is that the states are no longer describable as 
states with well-defined angular momentum, and the calculation of electron-phonon coupling is not straightforward. 
Therefore, we employ the EMA analysis of Takagahara for this. The Franck-Condon factors are estimated to be ~ 0.9 
between electronic states derived from the same multiplet. Calculations and experiments on Si reveal dark states with 
recombination times of microseconds fos!] . Tight binding states lack well defined quantum numbers. However, for 
spherical dots of 20 A there are multiple dark states which satisfy our phonon emission criteria |9^] . States from these 
multiplets can be used to form our logic and auxiliary states. Calculated Raman transitions between these states have 
quantitatively similar values to those obtained for CdTe above. 

Given an assumed radiative recombination rate F = 10'^ p5[| , we perform an analysis similar to the one above for 
CdTe. In Figure ^ A^max is plotted as a function of Luf for a variety of densities A and a threshold of one error every 
10 gates (e = 0.1). Figure ^ implies that one could construct a quantum computer composed of 700 quantum dots 
if e = 0.1. In Figure the extremum value of A^max, Nc is plotted as a function of e for a range of A values. The 
results are also summarized in Table 1. One sees that, unlike CdTe, for Si there does now exist the possibility of 
building a quantum processor that possesses an appreciably lower error rate of ~1 error every thousand gates. Most 
encouragingly, it seems possible to construct a small quantum information processor (5-10 qubits) with a larger linear 
support density IOAq, and an error rate of e < 10~^. Naturally, from an experimental perspective it would probably 
be more realistic to use a support having a density at least ten times greater than our proposed minimal density Aq 
that was estimated for a pure carbon chain (e.g., DNA Ea], carbon nanotubes, etched supports, etc.). 



VII. CONCLUSIONS 



We have developed a condensed phase scheme for a quantum computer that is analogous to the gas phase ion trap 
proposal and have explored the feasibility of implementing this scheme with semiconductor quantum dots coupled by 
a quantum linear support consisting of a string or rod. We have found that the Cirac-ZoUer scheme of qubits coupled 
by a quantum phonon information bus is also applicable in the solid state, and that there exist some advantages to 
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a condensed phase implementation. One such advantage is that there is a potential for significantly less noise in the 
information bus than in the corresponding gas phase scheme. Calculations by Roukes and co-workers ]10S| | suggest 
that much higher Q factors may be found for nanorods than are currently obtainable in ion traps. Clearly the extent 
of the usefulness of our proposal will be very dependent on the choice of materials. To that end we have analyzed the 
fidelity for two-qubit operations for several candidate systems, including both direct and indirect gap semiconductor 
quantum dots. We have presented the results of numerical calculations for implementation of the scheme with CdTe 
and Si quantum dots, coupled via either quantum strings or rods. While neither of these prototypical direct and 
indirect band gap materials reaches the level of fidelity and size required for large scale quantum computation, the 
indirect gap quantum dots (Si) do show a reasonably high fidelity with an array of a few tens of dots. 

One very revealing result of these explicit calculations of fidelity for one- and two-qubit gates is the limited scalability. 
The scheme initially appears highly scalable in concept due to the solid-state based architecture. However the detailed 
analysis given here showed that the dependence of the Lamb-Dicke parameter r] on the mass of the support is a basic 
problem which essentially limits the scalability to a few tens of qubits even in the more favorable indirect gap materials. 
The main drawback of this condensed phase scheme over the ion trap scheme is therefore the large reduction in rj 
deriving from the introduction of massive supports. Such a reduction has two important consequences. First, the 
laser intensities need to be increasingly large to perform operations faster than the decoherence time. Second, such 
large laser intensities necessitate the use of nodal and antinodal lasers |^,109|. Without these features, the probability 
of gate error is extremely high due to tran sitions to the carrier. This means that several of the alternative schemes 
proposed for ion trap computation ||lC , IIC | would not provide feasible condensed phase analogs (although the recent 
scheme of Childs and Chuang fo^ which allows computation with two-level ions (or quantum dots) by using both the 
blue and red sidebands could also be feasible in the condensed phase). 

Additional sources of decoherence which have been neglected here (scattering off the support, vibrational and 
electronic transitions in the support) will also act to limit the number of operations. However one source of decoherence 
which can be eliminated or at least reduced, is dephasing from the coupling to phonon modes of the support. This is 
a consequence of the requirement of extremely narrow band- width lasers, and therefore implies that a similar lack of 
dephasing will hold for other optical experiments on quantum dots which use narrow band-widths. One such example 
is the proposal to couple quantum dots via whispering gallery modes of glass microspheres p3| . More generally, this 
result offers a route to avoid dephasing for other spectral measurements on quantum dots ]114[ . 

An interesting additional application for this proposal is the laser cooling of nanorods. A single QD could be placed 
or even etched on a nanostructure. A laser tuned to the red support phonon side band of a QD excited electronic 
state would excite the energy of the nanocrystal, and at the same time lower the average phonon occupation of the 
support. When the unstable state relaxes, the most probable transition is the carrier transition. The net result is that 
the emitted phonon is blue shifted compared to the excitation pulse. The extra energy carried away by the emitted 
photon is thereby removed from the motional energy of the QD. 

The essential physical problem encountered in this condensed phase realization of the qubits coupled by phonon 
modes is the recombination lifetime of the qubit states, i.e., the exciton radiative lifetime. In principle this could be 
ameliorated by using hyperfine states of a doped nanocrystal. Recent experimental results dem onst rating electronic 



doping of semiconductor quantum dots offer a potential route to controlled access of these states [115|. The feasibility 



study presented in this paper does indicate that although the detailed physics of the qubits and their coupling is 
considerably more complicated in the condensed phase than in the gas phase, limited quantum computation may be 
possible with phonon-coupled solid state qubits. Further analysis and development of suitable nanoscale architectures 
and materials is therefore warranted. 
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IX. APPENDIX: COORDINATE REPRESENTATION OF ELECTRON AND HOLE STATES 

We give here the coordinate representation of the electron and hole states. These states were derived in (We 
employ a slightly different phase convention. Efros et al. when calculating the exchange Hamiltonian, ^Sh ■ Se, 
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between the hole and electron spin, Sh and S'e, use the convention, S,, ■ Se = S^S^ + ^{SJi;S- - Sj;S+). We instead 
use the convention that Sh - Se = S^S^ + ^{S'^S~ + S'^S'+) ). For convenience, we repeat the definitions of our qubit 
states [Eq. (^] in slightly more detailed notation: 

|0> = l*f/2,l/2(re))|*f/2,3/2(r/0> 
|l) = l^f/2,-l/2(re))|^3V3/2(r/.)> 

|2) = ^ [l*f/2,-l/2(l-e))|*f/2a/2(r/0> - l*f/2,l/2(l"e)>l*f/2,-l/2(l-h)) 

The electron states are simply the solutions to a free spin- 1/2 particle in a spherical hard- wall box: 



(46) 



where R is the radius of the dot, and Yj™ are spherical harmonics. 5* is a conduction band Bloch function and ±i is 
the z -projection of the electron spin: 



|5,±i) = |L = 0,mi = 0)|m, = ±i) 



(47) 



The holes states can be written explicitly as: 



l«'f/2,±l/2(r)> = -i?o(r)ro"l"±l/2) -i?2(r) (^y|5"2*>T3/2) + yf>^2^V±3/2) - 

= -i?o(r)roV±3/2) - i?2(r) (^y|r±>^i/2) - y|r2*V±i/2> - 

where Ri are the envelope functions and |wmj) are the valence band Bloch functions. 
The radial functions are: 



Y2>±l/2)j l*f/2,±3/2(r)) 



R2{r) 
Ro{r) 



A 
R^ 
A 



j,i^r/R) + -^^M^r/R) 



Joi^r/R) - 



where ji are spherical Bessel functions, [3 = mui/mfj.h is the ratio of the light to heavy hole masses, and is the first 
root of the equation 

JO(0)J2(V^0) + J2(0)jo(y^0) - 0. 

The constant A is defined by the normalization condition 

[Rl{r)+Rlir)] r^dr = 1. 

The valence band Bloch functions are given by: 

l"±3/2) = \L = l,mL = ±l)\ms = ±1/2) 
1 



l"±l/2) 



%/3 



V2|L = 1, mi = 0) |to, = ±1/2) + |L = 1, mi = ±1) |m, = Tl/2) . 

For our Raman transitions scheme, we have used states which were not analyzed in Jsot. I n particular, these states 
are from the 15elP5/2 exciton multiplct. To find these we used techniques developed in [116|, and then calculated the 
eigenstates of the exchange coupling using the method of [ |60[ . A Raman transition connects states of equal parity 
through a state of opposite parity. Therefore, the states of interest to us are the — ±1 state: 



18 



I*±l) = [l*f/2,T|l/2(re))|*f/2,±3/2(r/.)) + \/2|*f/2,±l/2(re))|*r/2,±l/2(r/.)) 

The electron state is as above. The hole state can be written explicitly as 



l*5''/2,±3/2(r)) = Riir) (^y|r°|«±3/2) + y|Y-i^V±i/2)^ 

+Rs{r) (^3yir3°|«3±/2) - iy|r3±>±i/2) + Y^n^VTi/2) + l^Y,^'\u^,^,) 



l*r/2,±i/2(r)) = Ri{r) (y^>'i^V±3/2)y|nV±i/2) + -/^n^^|«Ti/2)j 

where Ri are the envelope functions and \umj) are the valence band Bloch functions given above. 
The radial functions are 



^3(r) 
Mr) 



B 
B 



h{<t>'r/R) 



where cj)' is the first root of the equation 



and B is defined by the normalization condition 



/ {Rl{r)+Rl{r))r'dr = l. 
Jo 



Table 1. For a given error threshold, e, and support density A, the table shows the optimal value of -A^max,-^< 
CdTe and Si nanocrystals. Aq = 10 mu/A. 
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Laser 




Quantum 
Dot 



Linear Support 

FIG. 1. Schematic visualization of N quantum dots attached to a Unear support composed of a nano-scale rod or molecular 
string. Each quantum dot is addressed by a different laser. The absorption of the dots can be tuned by varying their sizes, 
allowing selective addressibility with lasers of different wavelengths. 
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FIG. 2. The addition of a sparse number of quantum dots to the linear support has relatively little effect on the dot modal 
displacement S„m.- Here we present results for Sn, the dot displacement resulting from the first harmonic of the support, 
for a system with two QDs attached to a string of length L = 2000 nm. The QDs are centered at 499 and 1501 nm. Each 
QD is represented as an increased density that is distributed over a length of 2 nm, e.g., 499 ±2 nm. Each dot experiences a 
displacement which is affected by the addition of the second dot on the support. The dot displacement measured relative to 
the value obtained from a homogeneous string (So), is plotted as a function of log((Ad + As)/As), where Ad is the linear density 
increment due to the dot and As the linear density of the string. The solid line guides the eye through the exact solution, the 
dotted line through the solutions obtained for a homogenous density equal to the average of Ad and As . One sees that the value 
of ySii changes by less than a factor of 2 over three orders of magnitude change in the normalized average density (Ad + As) /As. 
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FIG. 3. Energy level scheme for a quantum dot showing the laser fields and transitions necessary for one-qubit operations. 

(Energy level spacings are not to scale.) Levels |0) and |1) constitute the qubit. The auxiliary level |2) is not involved in these 
transitions and is not shown. The linear support modes are not involved cither. Two antinodal lasers, ki and k^, allow us 
to perform a Raman transition via a virtual state \v}. Transitions occur without changing internal phonon number, since the 
lasers frequency widths are smaller than the internal phonon frequency, ujf. 
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FIG. 4. Energy level scheme for a quantum dot on the linear support showing laser fields necessary for implementation of 
two-qubit operations. As described in the text, the use of nodal and antinodal lasers allows us to selectively transfer population 
from |0) to the lowest energy phonon sideband of |1) (labelled |l)|6i)) via a Raman transition, without transferring population 
to the carrier. The minimum phonon frequency is denoted cu'l. Non-resonant transitions to higher energy phonon sidebands, 
(|l)|s2)) constitute the main source of error in the proposed gates. For very high laser intensities, non-resonant quadrupolar 
transitions to higher level states (represented by the state \q}) also become important. 
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FIG. 5. Summary of energy scales involved in our proposal for single- and two-qubit operations. The states |0),|1), and |2) 
appear degenerate on the scale of this figure. 
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FIG. 6. Dependence of the maximum number of 20 ACdTe nanocrystal quantum dot qubits for which quantum computation 
is sustainable, subject to the three conditions determined by analysis of the two-qubit gate (see text): i) the fidelity per gate, 
T = .Fmax, ii) the antinodal laser intensity, I2 < I 2^^^, and iii) .7-"max > 1 — e, where e can be thought of as the error rate per 
gate frequency. The figure shows a plot of A'^max as a function of the frequency of the linear support phonon mode, ul, for two 
values of the linear support densities, Ao — 10 A = lOAo and e — 0.1. The extremum of the functions corresponds to the 

maximum possible scalability achievable for 20 A CdTe nanocrystal qubits. For larger values of ujI, the larger values of T^""^ 
which are possible in principle are offset by the need for higher intensity lasers. In this situation it is possible to support more 
qubits than are shown here by relaxing the first constraint. However one thereby loses the advantage of the increase in jTmax 
as is increased. 
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FIG. 7. Dependence of the optimal number of CdTe nanocrystals, Nc, (peaks in Fig ^ on the error threshold e plotted for 
various linear support densities. Aq = 10 




FIG. 8. Dependence of the maximum number of 20 A Si nanocrystal quantum dot qubits for which quantum computation 
is sustainable, subject to the three conditions determined by analysis of the two-qubit gate (see text): i) the fidelity per gate, 

= jFmax, ii) the antinodal laser intensity, I2 < I 2""^, and iii) J-"max > 1 — e, where e can be thought of as the error rate per 
gate frequency. The figure shows a plot of A^max as a function of the frequency of the linear support phonon mode, uf, for two 
values of the linear support densities, Xmin ~ 10 and A — lOAmin and e — 0.1. The extremum of the functions corresponds 
to the maximum possible scalability achievable for 20 A Si nanocrystal qubits. The degree of scalability is greater for the 
indirect band gap material than for the direct band gap CdTe nanocrystals shown in Figure ^, and shows less dependence on 
the linear support density A. 
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FIG. 9. Dependence of the optimal numer of Si nanocrystals, Nc, (peaks in FigM) on the error threshold e plotted for various 



linear support densities. Aq = 10 
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